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FOREWORD 

This document represents the final report of a study 
performed by Lockheed Missiles & Space Company, Inc., 
Huntsville Research & Engineering Center, entitled "Hyper - 
velocity Meteoroid Impact on Orbital Space Stations, " Con- 
tract NAS8-28473. The program was administered by the 
National Aeronautics and Space Administration's George C. 
Marshall Space Flight Center under the direction of Mr. 

David W. Jex, Space Sciences Laboratory. 

The authors are grateful for the guidance and assistance 
provided by Mr. Jex in addition to Dr. L. Raman, Mr. R. L. 
Holland and Mr. G. Heintze. 
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SUMMARY 

In the event of a hypervelocity impact of a meteorite on a spacecraft, 
structural damage may result. Of particular interest is the backside spalla- 
tion caused by such a collision. To treat this phenomenon two numerical 
schemes were developed in the course of this study to compute the elastic- 
plastic flow and fracture of a solid. The numerical schemes are a five-point 
finite difference scheme and a four-node finite element scheme. The four- 
node finite element scheme proved to be less sensitive to the type of boundary 
conditions and loadings. Although further development work is needed to 
improve the program versatility (generalization of the network topology, 
secondary storage for large systems, improving of the coding to reduce the 
run time, etc.), the basic framework is provided for a utilitarian computer 
program which may be used in a wide variety of situations. Analytic results 
showing the program output are given for several test cases. 
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NOMENCLATURE 


Symbol 

A 

a l» a 2’ a 3* b l* b 2 



m 


P 

Q 

2 

q 

q' 

s 

t 

U, V 

V 


Description 

surface area 

coefficients of equation of state 

damping coefficient for quadratic damping 
scalar damping coefficient 
energy of the control volume 
internal energy per unit mass 
any function 

body force, body force per unit mass 

surface force 

nodal forces 

mass 

pressure 

heat transfer rate per unit area per unit time 

heat transferred into control volume 

viscous deviator stress tensor 

viscous pressure 

position vector 

time 

displacements 

volume 


v 
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Symbol 

W 

x, r, 0 


Greek 

T 


6 



6 

rx 


€ 


r) 

X 

M 

Mi 


P 

CT 

a r a 2 ,a 3 

2 

U) 


NOMENCLATURE (Continued) 

Description 

work done on control volume 
cylindrical coordinate system coordinates 
yield modulus 


unit dyad (tensor) 

coordinate system rotation corrections 
strain 

density/rest density ratio 

plastic flow correction factor 

shear modulus 

viscosity coefficient 

0 is planar motion, 1 is axisymmetric 

density 

distortional stress deviator tensor 
principal stress 

stress tensor 
rotation angle 


Superscripts 

• 

e 

n 


denotes time differentiation 
designates element 
denotes time step 


Subscripts 

x, r, 0 

xr 

o 


denotes component in coordinate direction 
off diagonal term in stress and strain tensors 
denotes rest condition 
vi 
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Section 1 
INTRODUCTION 

The analysis of dynamic response of elastic media to rapid loading, in 
particular meteorite impact, has been of considerable interest to the scientific 
community. Some of the most notable analytical investigations are given by 
MacCormack (Ref. 1), Rosenblatt (Ref. 2); Read and Bjork (Ref. 3) and Wilkins 
(Ref. 4). In these studies the emphasis has been placed on the dynamic response 
of the target or impacted material and on the response of the projectile. 

Experimental studies (Ref. 5) have shown that for a given weight, a thin 
meteorite bumper arrangement gives better protection than a single skin. The 
purpose of the bumper is not necessarily to completely stop the incoming pro- 
jectile but to at least absorb some of the incident energy before failing under 
the applied load. Backside spallation of the bumper then scatters fragments 
of both the meteorite and the bumper material which in turn strike the structural 
skin. These fragments possess a lower total energy and also impact the struc- 
tural surface over a much wider area thus decreasing the peak loading. 

Development of the bumper approach created a need for an analysis 
which treats the fracture phenomena of a material. It is the purpose of this 
study to provide this capability. The schemes chosen are based on combina- 
tions of ideas found in the above references and the work of Ang (Refs. 6, 7, 8). 

Of the two schemes investigated, a five-node finite difference scheme and a 
finite element scheme, the finite element scheme produced the best numerical 
behavior. It can be shown that the finite element treatment actually is analogous 
to a nine-point finite difference scheme. The additional four points enhance the 
stability characteristics (with respect to the five-point scheme). 

In the following discussion both approaches are discussed. Computer 
programs have been written to conduct the calculation and are briefly discussed 
in this report. 
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Section 2 

EQUATIONS OF MOTION GOVERNING 
CONTINUOUS MEDIA 


For a control mass in motion conservation of mass becomes 


dm _ d C 

dt dt J 1 

cv 


pdV = 0 


from Newton’s second law we get 


~ (mS ) = f dF + f dF 
dt me J s I B 

c s cv 


while the first law of thermodynamics yields 


dE dD dW 


dt dt dt 


Now the velocity of the mass center is 


me 


= J' pS dv/m 


so that Eq. (2.3) becomes 


cv 


-£_/ pSdV J dF s *f dF 


B 


cv 


cs cv 


The work done on the control volume in time dt is 


AW =J dF g • dS +y* pf B • dS dV 


c s 


cv 
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where dS is the distance moved in that time so that the instantaneous rate of 
work done on the control volume is 


W = J S • dF g + / pS • 




cs 


cv 


The rate of change of energy contained within the control volume is 


E =/p!e + i S. S JdV 


CV 


It is convenient to express 


d 2 f 

ar “/ 


Q ’ dA where Q is the heat transfer 


cs 


vector per unit time per unit area. In a similar fashion let 


dF = 2 • dA 
s 


so that the governing equations can be summarized as below 


4 / 


pdV = 0 


cv 


(2.4a) 


&/pSdV =/ S .dA + J pf B dV 


CV 


cs 


cv 


(2.4b) 


dl ft 


cv 


p (e + \ S • S) dV =y* (Q + S • 2) • dA + / pS • f R dV 


CS 


J P i 


B 


cv 


(2.4c) 


where the control surface is determined by the space time history of points on 
the original or starting surface. 
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Assuming continuous functions within the region we may write, with the 
aid of Gauss' theorem, for momentum and energy 

pS = V * S + pfg 

p (e + S • t) = V • (Q + S • S) + pt • I B 
Combining (2.5a) and (2.5b) yields 

• 

pe = V • Q + (S • V) • S (2.6) 

At this point it is convenient to decompose the dyad (or second order tensor) Z 
into viscous and inviscid components viz: 

S=Sj+S v = a - pd - q'^5 - q 

Now a is the deviator stress tensor due to elastic/plastic distortions, p is the 
hydrostatic pressure, q' the viscous pressure and q is the deviator viscous 
stress tensor. The elastic stress convention is such that tension yields a posi- 
tive stress while pressure in hydrodynamics is usually defined in the opposite 
sense. The viscous stresses usually are defined in the same sense as pressure. 


(2.5a) 

(2.5b) 


Assuming an axially symmetric motion (or a planar one) the equations of 
motion may be written, 


M 

pr 



as 


rx 

ax 



px 



(2.7a) 


(2.7b) 


for a cylindrical coordinate system. 
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The equation of continuity becomes 


V 

V 


e + e 
r x 


+ Ve, 


(2.7c) 


where the strain rates are; 


* - 9 r _ 9 v 

€ r 9 r 9 r 

• _ 9 x _ 9 u 

€ x 9 x 9 x 


( 2 . 8 ) 


r _ v 
r r 


€ 

rx 


9 r , 9 x 

9 x 9 r 


9 v , 9 u 
9 x 9 r 


The material law for linear elastic isotropic homogeneous material is 

a =2 ix (e 


I V. • 
r 3 V ^ r 


, IV.- 

" * € x " 3 V * + 


<j Q = Zy. (e 0 - (axisymmetric only) 


(2-9) 


cr = a e +6 
rx ^ rx rx 


The quantities 6,6,6 are corrections to account for the rotation of the 
r x rx 

coordinate system. 

Consider the rotation of the coordinate system by an angle u. Then 


2 

6 = (2 - 2 ) sin u - 22 sinu cosco 

r ' x r' rx 

6 = -6 


( 2 . 10 ) 


rx 


-22 sin u) - (2 - 2 ) sinw cosu 

rx x r 
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but 


w = V x (ui + vj) 


Equations (2.9) could be differentiated with respect to time but there is no need 
since the integral is the quantity that is really required. 


Expansion of Eq. (2.6) yields 


pe 


2 e + 2 
xx r 


€ + 
r 


e e 


+ S'' 
xy 


e 

xy 


which yields 


pe = -(p+q’) ^ + ( a x " € x + (a r - q r ) e r + (a ( 


-V 6 0 


( 2 . 11 ) 


+ (ct - q ) e 
' xr ^xr xr 


The hydrostatic pressure must be obtained from an equation of state. 
Functionally we may express this as 

p = p(p, e) 

The form used by Wilkins (Ref. 4) is 

p = a^ (rj - 1 ) + a^ (*] - 1 )^ + 83(11 - 1 )^ + £bj (n- 1 ) + b ^ (*? - 1 )^J e 

where 11 = p/p Q and is the functional used in this study. 

To complete the governing equations the viscosity must be discussed. 
Generally speaking in numerical analysis the introduction of viscosity occurs 
more frequently to enhance stability of the numerical solution than for the 
purpose of treating internal damping or viscous effects. As such the 
viscosity models are created to control some undesirable behavior and do not 
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necessarily have firm underlying mathematical principles. Several mathe- 
matical models were tested in the subject analysis. These are: (1) Linear 

Damping; (2) Quadratic Damping; (3) Velocity Angular Distortion; (4) Spatial 
Angular Distortion; and (5) Spatial Displacement Distortion. 

The first method has firm mathematical grounding in the sense that it 
is a conventional model for fluids obeying the Reiner- Revlin law. The second 
method uses a first coefficient of viscosity which is proportional to -the diver- 
gence and a second coefficient which is zero. The last three methods are 
artificially contrived. The third was suggested by Wilkins (Ref. 4) to control 
an undesirable "hourglass" distortion. The last two were developed in this 
study due to the inadequacy of the five point lattice scheme to "remember" 
what the original rest configuration was and to return to that state in the 
absence of disturbing forces. 


In the scalar damping 

• ^ V 

q' = c gMl y 

q r = 2m 1 * € r ‘ I V * 
= 1 * C x ‘ 3 V ^ 

q e 2/i i * € e " 3 v ) 

q rx ^ 1 ^ rx 


(2.12a) 
(2.12b) 
(2.12c) 
(2. 12d) 
(2. 12e) 


while in the quadratic damping 


q'. 


2 V 2 

C oPo<V> 


(2.13) 


A departure from Wilkins is taken in that the factor was omitted in the 
quadratic damping since it precluded a uniform radial distribution in cases where 
no variations of properties in the radial direction can be expected. 
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The remaining topics of yield and fracture complete the governing 
equations discussion. The principal stresses are 


and 


1.2 


a + a 
r x 


2 







+ a ) 

x 


(2a )‘ 

rx 


(2.14) 


a^ = a^ (axisymmetric only) 

The yield condition is given by Von Mises as 

2 . 2,2 2 . . 2 . 

a l + CT 2 + CT 3 ' 3 (Y o ) - 0 

If this condition is not satisfied then all stresses are corrected by multiplying 
by the factor 



The validity of this procedure is discussed by Hill (Ref. 9). A fracture occurs 
if the principal tensile stress exceeds the ultimate tensile strength. 

This concludes the discussion of the governing equations for the continuous 
media. The numerical analogs are developed in the following section. 
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Section 3 

NUMERICAL ANALOGS TO THE GOVERNING EQUATIONS 

3.1 FIVE-NODE FINITE DIFFERENCE SCHEME 

This is the simplest numerical analog possible, at least with regard to 
the region of influence. For each node point in the array the positions, dis- 
placements, stresses, pressure, viscosity, and vel ocities are stored. Figure 
1 illustrates a typical interior node. 

4 

1 

3 

Fig. 1 - Typical Five-Node Lattice Numbering Scheme 
The spatial derivatives of any function f are 

jfx = 2A | * f 2 V * r 4 ' r 3^ ' ^ f 4 " V * r 2 ” r l^| (3.1.1) 

8~r ‘ 2A { ^ f 2 " f l^ ^ X 4 ' x 3^ ‘ ^ f 4 " f 3^ ^ X 2 ' x l^ j (3.1-2) 

The accelerations are computed using Eqs. (2.7a) and (2.7b). Integration of 
the accelerations yields new velocities over the entire field while a second 
integration yields the new positions and displacements. Using the above 
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information the strain rates at the new time may be found from Eqs. (2.8) while 
the new deviator stress distributions are found by integrating Eqs. (2.9). 

The pressure and energy may then be evaluated and the new viscous terms 
determined. From this information the updated stress tensor is found. 

A special matrix is retained which for each node point stores the index 
of the point to the left, right, top and bottom. In the event that point 0 is on 
the left boundary the index for point 1 is set to that of point 0 and a one sided 
derivative results. This process is used for all peripheral points. As will 
be seen later the unconstrained corner point such as shown (Fig. 2) has poor 
stability and constitutes the greatest problem with the five-point method. 


+2 


* 3 

Fig. 2 - Degeneration of Five-Point Grid at Corner Point 


Boundary conditions are treated with special subroutines which reset displace- 
ments, velocities and stresses in an appropriate fashion. These routines are 
extremely flexible such that any type of boundary conditions may be treated by 
recoding rather than by data card input. Very few cards are involved, typically 
5 to 10, in this recoding. This approach was considered superior to a complex 
input setup in which every conceivable type of boundary condition was anticipated. 

3.2 FOUR -NODE FINITE ELEMENT SCHEME 

This computer program features some of the aspects of the traditional 
finite element programs used in structural analysis. Dynamic core storage 
allocation is used for the major size -dependent variables, which are grouped 
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in nodal quantities (position coordinates, displacements, velocities, forces, 
masses) and regional quantities (stresses, pressure, scalar viscosity, total 
energy and distortional energy). 

The elastic -plastic continuum is replaced by four-node regions (elements) 
as shown on Fig. 3. The corner points are numbered 1 through 4 counterclock- 
wise. The state of stress in the element is assumed to be constant. „ 

The mass characteristics of the continuum are modeled as discrete lumps 
located at the node points (Fig. 4). The motion of the continuum is treated in 
terms of the motion of these lumps. 

The equation of motion is simply a sum of all forces, elastic and inertia, 
acting on the lump. No special logic is required for boundary lumps since 
the elastic forces are summed as the element stresses for each element are 
updated. This simplifies the coding considerably, doing away with special 
features for various types of boundaries. The equation of motion may be stated 
as 


mx = ^ F® (3.2.1) 

e 

mr =2 F r (3.2.2) 

e 

e e 

where F x and F^ are the elastic (spring) forces of the elements holding lump 
m. 


The velocities are obtained by taking 



(3.2.3) 
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^ A, X 


Typical Node With Surrounding Elements 
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S f; a«" +1 


(3.2.4) 


If a velocity boundary condition is imposed the affected nodal velocities com- 
puted from Eqs. (3.2.5) and (3.2.6) are simply reset to the imposed values. Force 
boundary conditions, such as externally applied forces, are imposed before Eqs. 
(3.2.5) and (3.2.6) are executed. 


In another integration step the displacements are computed, 


n+1 n . «n+T A .n+T 

u = u 4 x ^ At ^ 


n+1 n , -n+i A n+i 

v =v+r ^ At z 


Displacement boundary conditions may also be imposed. 


(3.2.5) 

(3.2.6) 


Now the element quantities can be computed based on the new velocities 
and displacements, from the equations of motion. 


A fracture condition has been added, which severs the springs representing 
the strength of the element. This fracture occurs when the principal tensile 
stress exceeds the ultimate tensile strength of the material. 


The corner forces are computed for the element based on the new stresses 
and damping terms. The directional damping terms are not stored. These 
corner forces are 


F xl * -^[' 2 xl (r 2- r 4>- 2 xr ,X 2- !t 4>] 1:+ ? S xr A | t 
F x 2 = -i [< 2 x2( r 3- r l)- 2 xr (lt 3- x l»] F+ ? S xr A|,! 
F x 3 * -i[( S x3 (r 4 ' r 2> - S xr (x 4 - x 2>] 7 + 5 Z xr A l* 
F x4 = -i [< Z x4 ('1 - ’• 3 ' - 5cr < X 1 - x 3>] ? + Kr A ^ 
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F rl “ i [ Z yl < x 2 - x 4> - Z xy < r 2 ‘ r 4 >J r + 5 < Z r ‘ Z 8> A I* 

F r2 = *[ Z y2 (x 3 ' x l* " Z xy (r 3 ‘ r l*] r + ? (S r " Z e ) A I* 

i r l_l I * (3.2.10) 

F r3 = * [ Z y3 ,X 4 - x 2> - Z xy < r 4 ’ r 2>J r + 4 <2 r ‘ Z 9 > A 1 

F r4 = *[ Z y 4 < X 1 - x 3> - Z xy < r l - r 3>] 7 + ? <=r ‘ Z 0 > A I* 

I* 

r = 1; I =0 for plane strain problems. 
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Section 4 

NUMERICAL RESULTS 

To investigate the numerical problems involved in the computation of 
elastic -plastic flow a one -dimensional code was written and implemented. A 
wave propagation problem given by Wilkins (Ref. 4) was checked with this code 
using two different mesh sizes (250 and 50 node points). The material is de- 
scribed by P = 0, p^ = 4g/cm p = 3 Mb and Y q = oo. The left boundary of the 
^5 cm thick plate is excited by x = -10"^ sin2?rt cm/p sec. 

The stress waves are depicted on Figs. 4-1 through 4-3. 

The five -node finite difference scheme was first tested on a one -dimen- 
sional problem although the calculations were performed two-dimensionally. 
Quadratic damping was used in this test case. A uniform velocity was imposed 
at one end of a rod while the other end was fixed. The rod was constrained to 
slide in a tube such that no vertical displacement was possible. The axial 
velocities resulting from the improved boundary conditions are shown in Fig. 
4-4 for several node points. As the time approached infinity the internal 
oscillations damped and a linear velocity distribution throughout the rod was 
achieved. 

The two-dimensional codes initially suffered from stability problems. 

The examples shown were taken again from Wilkins, with the material data 
P = 1.88 (il - 1) Mb, p Q = 7.72 g/ cm^, p = 0.814 Mb and Y = oo . The cantilever 
plate is 1.00 by 5.25 cm, the fixed-fixed plate is 1.00 by 10.50 cm. The excita- 
tion was achieved by supplying a boundary velocity. 

Figures 4-5 and 4-6 show how the upper edge of the cantilever plate is 
set into motion at constant velocity and constant rotation, respectively. Rela- 
tively smooth deformation patterns result. If an edge is released, however, 
instability occurs almost immediately (Fig. 4-7). This points to a weakness 
in the formation in the edge condition of the five -point scheme. For all three 
plots the deformations are magnified ten-fold. 
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Fig. 4-1 - Stress Distribution in the Plate 
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Axial Distance (cm) 


Fig. 4-2 


- Reflection of the Wave at the Free Boundary 
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Fig. 4-3 - Stress Distribution for Coarse Mesh at t = 3.3 n sec 
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VIBRATION OF AN ELASTIC PLATE CLAHPEO AT ONE END 



Fig. 4-5 - Elastic Shear Deformation (Five-Point Scheme, Plane Strain) 
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VI0BATION OF AN OLA S1IC PLATE CLAMPEO AT ONE ENO 



TIHf MW mCROSEC ' CVCLE NUMBER 

Fig. 4-6 - Elastic Bending Deformation (Five-Point 
Scheme, Plane Strain) 
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VIBRATION OF an ELASTIC PLATE CLAWED AT ONE END 



Fig. 4-7 - Instability of the Five-Point Scheme at 
the Corner Points 
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The fixed -fixed plate problem was solved by the five -point scheme, but 
shear waves occurred (Fig. 4-8). These could neither be controlled by the 
angle q nor by the Navier-Stokes q. Several corrective schemes were tried. 
The stabilizing effect was good for considerable length of time but eventually 
instability occurred. 

For rotationally symmetric problems the five-point scheme is consider- 
ably more stable. Figures 4-9 and 4-10 show a circular plate with a hole 
clamped at the outer edge. It is significant that the corners on the inner 
edge remain stable. Eventually the integration becomes unstable, however, 
although the grid remains smooth. 

Results obtained with the four-node element scheme have shown good 
comparison with those reported by Wilkins as far as amplitude and frequency 
of the vibrating plate is concerned (Fig. 4-11). Unfortunately an hourglass 
distortion appears. This should be controllable by a higher coefficient for the 
angle damping coefficient. 
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VIBRATION OF AN ELASTIC PLATE CLAMPED AT BOTH ENOS 



TTME» tt.683 HICROSEC CYCLE NurBER 120 


Fig. 4-8 - Shear Waves in Elastic Plate (Five-Point 
Scheme, Plane Strain) 
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VIBRATION or AN ELASTIC PLATE CLAHPEO AT CIRCUMFERENCE 



Fig. 9 - Cutoff at t = 8ft sec 
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VIBRATION OF AN ELASTIC PLATE CLAMPED AT CIRCUMFERENCE 
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Fig. 4-10 - The Grid Moves to the Other Side 
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VIBRATION Of A CANTILEVER PLATE 



TINI« ITTjlll NICROBEC CTCLE NUNBER TIB 


Fig. 4-11 - Bending of a Cantilever Plate Showing 

Hourglass Distortion 
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Section 5 

CONCLUSIONS AND RECOMMENDATIONS 

In order to analyze the behavior of a spacecraft structure subjected to 
a meteorite impact two numerical analogs to the governing differential equa- 
tions related to elastic-plastic flow and fracture of a material have been 
developed and computer programs have been written to perform the calcula- 
tions. Of these two schemes, the finite element or four-node scheme appears 
to be superior for the following reasons. 

• Stability problems appear to be only of the "hourglass" 
type and may be controlled. 

• It is simpler to formulate 

• The boundary condition treatment is simple to apply. 

Full production versions have not been created with either method. Stability 
problems occupied more effort than originally anticipated and thus detracted 
from development of a user- oriented analysis. 

It can be said, however, at least regarding the finite element scheme 
that the major problems have been overcome and that a nominal amount of future 
development will, in fact, produce a utilitarian program. 

During the course of the study the fracture condition itself was not 
thoroughly investigated, but no undue problems are anticipated. 

It is concluded that further exercise of the existing programs and further 
development is warranted. 
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Appendix 

USERS' MANUALS FOR THE COMPUTER PROGRAMS 
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A. 1 FIVE -POINT SCHEME 


Card 

Columns 

Format 

Description 

1 

1-72 

12A6 

Heading 

2 

1-5 

15 

Time step increments for output 


6-10 

15 

Total number of time steps 


11-15 

15 

Type of symmetry 0 = plane 

1 = cylindrical 


16-20 

15 

Number of stations in x-direction 


21-25 

15 

Number of stations in r -direction 


26-30 

15 

Type of output 0 = plots desired 

3 

1-10 

E10.4 

p (g/crn ) = density 


11-20 

E10.4 

p (Mb) = shear modulus 


21-30 

E10.4 

Y° (Mb) = yield strength 


31-40 

E10.4 

(Mb) = ultimate tensile strength 


41-50 

E10.4 

C 2 = a constant 
o 

4 

1-10 

El 0.4 

Coefficients of equation of state 


11-20 

E10.4 

Coefficients of equation of state a 2 


21-30 

El 0.4 

Coefficients of equation of state aj 

5 

1-10 

El 0.4 

fy = length of the region in r-direction 


11-20 

E10.4 

fx = length of the region in the x-direction 


21-30 

E10.4 

y a = beginning of region in r-direction 
(plate with hole) 


31-40 

E10.4 

a = magnification factor for plots 

6 

1-72 

12A6 

Message to the operator of the SC 4020 
plotter. 

A. 2 

FOUR-NODE ELEMENT SCHEME 


Card 

Columns 

Format 

Description 

1 

1-72 

12A6 

Heading 

2 

1-5 

15 

Time step increments for output 


6-10 

15 

Total number of time steps 


11-15 

15 

Type of symmetry 0 = plane 

1 = cylindrical 


16-20 

15 

Number of stations in x-direction 


21-25 

15 

Number of stations in r-direction 


26-30 

15 

Type of output 0 = plots desired 
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Card Columns Format Description 

3 1-10 E10.4 p (g/cm^) = density 

11-20 E10.4 n (Mb) = shear modulus 

21-30 E10.4 F w (Mb) = yield strength 

31-40 E10.4 F^ (Mb) = ultimate strength 

41-50 El 0.4 C Q ' 

51-60 E10.4 C = constants 

61-70 E10.4 C^ J 

71-80 E10.4 a g (cm/p sec) = speed of sound 

4 1-50 5E10.4 Coefficients of equation of state 


5 

1-10 

E10.4 

X A 


11-20 

E10.4 

X B 


21-30 

E10.4 

y a 


31-40 

E10.4 

y b 


Boundaries of the grid 


6 1-72 12A6 Message to the operator of the SC 4020 

plotter. 
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